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Only a severely restricted class of tensor fields can provide classical spacetime geometries, 

^ ' namely those that can carry matter field equations that are predictive, interpretable and 

I , quantizable. These three conditions on matter translate into three corresponding algebraic 

5-1 . 

conditions on the underlying tensorial geometry, namely to be hyperbolic, time-orientable 
and energy-distinguishing. Lorentzian metrics, on which general relativity and the standard 
^ \ model of particle physics are built, present just the simplest tensorial spacetime geometry 

^\ , satisfying these conditions. The problem of finding gravitational dynamics — for the general 

■ tensorial spacetime geometries satisfying the above minimum requirements — is reformulated 
^SJ ■ in this paper as a system of linear partial differential equations, in the sense that their 

O! 

_ solutions yield the actions governing the corresponding spacetime geometry. Thus the search 

" . l , for modified gravitational dynamics is reduced to a clear mathematical task. 
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I. INTRODUCTION 

Over the past two decades, remarkable gaps in our understanding of matter have been revealed — 

n 

astrophysical observations [1] by now overwhelmingly indicate that only about 4% of the matter 
and energy in the universe can be of standard model origin. Indeed, one cannot attribute the 
remaining 21% of dark matter or 75% of dark energy to standard model matter or its vacuum 
energy. For on the one hand, does the bullet cluster 3] show the existence of only gravitationally 
interacting dark matter. On the other hand does the calculation of dark energy as the vacuum 
energy of standard model fields yield a result that is infamously off the observed value by 120 
order of magnitude Q], which is jokingly referred to as the worst prediction of elementary particle 
physics. It is consensus that these observations thus point at something fundamental we do not 
understand about matter or gravity. 

Going deeper than just postulating modified Lagrangians for either matter or gravity, however, 
one quickly realizes that modifying matter and gravity dynamics independently quickly becomes in- 
consistent. This is because they both build on — and more importantly, are both tightly constrained 
by — the common underlying spacetime geometry. 

For the Einstein equations, on the one hand, present the unique dynamics with a well-posed 
initial value problem one can give to a Lorentzian manifold [J]. Thus modifying gravitational 
dynamics necessarily comes at the cost of deviating from Lorentzian geometry as the spacetime 
structure. This is illustrated for instance by Brans-Dicke gravity and its avatars as well as f{R) 
gravity theories, which all feature at least an additional scalar gravitational degree of freedom. But 
then one needs to couple matter fields to the corresponding modified spacetime geometry, and one 
sees that modifying gravitational dynamics compels one to think about modified matter dynamics. 

Vice versa, even minimal deviations from the standard model of particle physics on the other 
hand quickly produce matter dynamics whose causality does not coincide with the causality defined 
by the underlying Lorentzian manifold. Famously, this surprisingly already happens for some fully 
covariant Lagrangians [5]. If one wishes to consider such matter dynamics on a Lorentzian manifold, 
let alone any even slightly more exotic matter, one requires also modified gravity dynamics to 
provide backgrounds that make the matter equations causal. Thus unless one entertains the claim 
that all matter that could exist in Nature must be made precisely in the image of Maxwell theory — 
which in the face of 96% of all matter and energy in the universe being of entirely unknown origin 
appears an unnecessarily restrictive and arbitrary idea — one sees that modifying matter dynamics 
generically compels one to construct a corresponding modified gravity theory to be consistent. 
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Therefore if one sets out to modify the otherwise time-tested theory of Einstein gravity coupled 
to standard model dynamics, the intimate link between consistent matter dynamics and gravi- 
tational field equations compels one to be particularly careful, and base any modification on a 
common underlying geometry that is so constrained as to make the combined theory work. At 
the very least, it is to be physically required that the spacetime geometry render the entirety of 
observed matter field dynamics coupling to it predictive, interpretable and quantizable [3]. The 
identification of all tensorial geometries complying with these minimum criteria and the determi- 
nation of their gravitational dynamics is the purpose of the present paper. 

Fortunately, these rather fundamental physical conditions translate into three simple algebraic 
conditions Q] that an otherwise arbitrary tensor field must satisfy in order to provide a valid 
spacetime geometry: it must be hyperbolic, time-orientable and energy-distinguishing, as we will 
explain in the first technical section below. Thus the spectrum of tensor fields that can serve as a 
spacetime structure in the presence of specific matter field dynamics [25] is greatly restricted. So 
restricted, in fact, that all kinematical constructions familiar from the special case of Lorentzian 
metrics can be made, for precisely the same structural reasons, also for any such tensorial spacetime. 

The deeper dynamical principles behind Einstein's field equations, which were revealed by 
Arnowitt, Deser and Misner 7^ a long time ago, are fortunately not swept away with a change 
from Lorentzian geometry to one of the alternative tensorial geometries described above. From 
the geometrodynamical point of view, gravitational dynamics is all about evolving the spatial 
geometry from one suitable initial data surface to an infinitesimally neighbouring one, such that 
ultimately all spatial geometries recombine to an admissible spacetime geometry; indeed, applying 
this principle, Hojmann, Kuchar and Teitelboim J, ^ derived the Einstein-Hilbert action with 
numerically undetermined gravitational and cosmological constants as the unique dynamics for the 
special case when the tensorial geometry is in fact a Lorentzian metric. But this geometrodynamic 
principle stands for any tensorial spacetime geometry 26|]. And it is the purpose of the present 



article to show that the very same principle indeed determines the dynamics of any hyperbolic, 
time-orientable and energy-distinguishing tensorial spacetime geometry. But indeed only of such. 

The main result of this article, beyond its technical details, is thus the observation that the 
search for gravitational dynamics beyond general relativity can be reduced to solving a mathemati- 
cal representation problem. This is achieved by invoking precisely the same geometrodynamic prin- 
ciples as followed by jS], but applied to any hyperbolic, time-orientable and energy-distinguishing 
tensorial spacetime geometry. As a consequence, gravitational dynamics (including and beyond 
{2?! Einstein) need not be postulated, since they can be derived for any tensorial spacetime geom- 
etry, as we will show, by solving a a system of homogeneous linear partial differential equations. 
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Thus the question whether there is an alternative to general relativity as a classical gravity theory 
carries over in the formalism developed in the paper to a mathematical existence problem. Now if 
such solutions exist at all, the question of whether there are other dynamics for the same geometry 
translates into the mathematical question of the uniqueness of such a solution. And finally, the 
problem of constructing concrete alternative gravity dynamics amounts to nothing less, but also 
nothing more, than finding an actual solution of the linear system of partial differential equations. 
In the philosophy of this paper, studying modified theories of gravity therefore amounts to studying 
solutions to the said equations. 

The organisation of this paper is as follows. In section Hi Al we start by giving a concise review 
of tensorial spacetime geometries, culminating in the insight of how normal co-vectors to initial 
data surfaces are mapped to normal vectors, which is the pivotal technique in constructing the 
canonical dynamics for such geometries. This is then used in section Hi B I to derive the deformation 
algebra of initial data surfaces in any tensorial spacetime geometry. Since geometrodynamics 
must evolve geometric initial data from one initial data hypersurface to the next, gravitational 
dynamics must represent this deformation algebra on a geometric phase space, which we explain 
in section III CI Chapters IIIII and IIVI then deal with different incarnations of the same program; 
the former derives the dynamics of dispersion relations independent of any underlying geometric 
tensor, while the latter deals with the more fundamental question of deriving dynamics for the 
fundamental geometric tensor. More specifically, the supermomentum for the dispersion relation 
geometry is constructed in section IIIIBl while the corresponding superhamiltonian splits into a 
non-local part that we construct in section IIIICI and a local part that is determined by equations 
that we derive in the course of sections IIIID| IIIIEI and IIIIFI Section IIIICI then derives further 
insight on the structure of the equations determining the local part of the superhamiltonian for 
dispersion relations, before we recover the dynamics for the dispersion relation of standard general 
relativity in section IlII HI starting from nothing more than the dispersion relation of lowest possible 
degree, demonstrating that the general principles underlying our study of all spacetimes are none 
other than those also underlying general relativity. The construction of the supermomentum and 
superhamiltonian determining the dynamics of a fundamental geometric tensor field in chapter 
IIVI proceeds very much along the same lines, but depends heavily on the algebraic structure of 
the fundamental geometric tensor one considers. Building on some work recycleable from section 
IIIII we derive the equations whose solution yields the dynamics for area metric manifolds as a 
prototypical example in section IIVCI after having constructed the relevant geometric phase space 
in section IIVBI We conclude in section |Vl 
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II. KINEMATICS OF TENSORIAL SPACETIMES 

Whether a tensor field can provide a spacetime structure depends on the matter one wishes to 
consider on it. In this chapter we will first review how three fundamental requirements one needs to 
ask of any realistic matter theory — predictivity, interpretability and quantizability — greatly restrict 
the tensor field backgrounds they can couple to. In particular, we will see that the dispersion 
relation associated with the entirety of field equations on a spacetime plays a prominent geometric 
role, and almost single-handedly encodes the kinematics associated with the underlying tensorial 
geometry. These insights, which are reviewed here in section III Al in a brief but self-contained 
manner for the benefit of the reader, are then used in section III Bl to study the deformation algebra 
of initial data hypersurfaces in tensorial spacetimes. The basic geometrodynamic idea, namely 
to use this deformation algebra to derive canonical dynamics, is then laid out in section III CI and 
presents the key to constructing gravitational dynamics for the dispersion relation in chapter IIIH 
or more fundamentally, the underlying tensorial geometry in chapter IIVI 



A. Primer on tensorial spacetime geometries 



In this section, we give a concise review of tensorial geometries that can serve as a spacetime 
structure. The technical proofs underlying this summary are presented in detail in ^ and rather 
pedagogical fashion in the lecture notes [9]. To aid the reader's intuition, we illustrate each abstract 
construction in this section immediately for the familiar example of a standard metric geometry, 
before mo ving on to the next construction. Occasionally we will also contrast this to area metric 



geometry 



IC 



111 ] as a comparatively well-studied example for a non-metric tensorial geometry. 
Having studied the general theory and these examples, the reader should be in a position to carry 
out a similar analysis for his or her favourite tensorial geometry. 

. So we consider, in 



All we know about spacetime, we know from probing it with matter 
addition to an a priori arbitrary tensor field G (the "geometry") on a smooth manifold M also 
a field (j) (the "matter"), which takes values in some tensor representation space V and whose 
gauge-fixed dynamics are encoded in linear field equations that transform as a tensor. Since the 
only other structure on the manifold besides the matter field is provided by the geometric tensor 
G, the coefficients Q of the matter field equations (after removing potential gauge symmetries and 
separating off the related constraint equations) must be built entirely from the geometric tensor 
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and its partial derivatives 29|], 

N 

J] Q[G]%-^- {X)da, . . . da„<P^{x) = , (1) 
n=0 

with small latin spacetime indices running from 0, . . . , dimM — 1 and capital latin respresentation 
space indices ranging over 1, . . . ,diml/. It is straightforward to establish that in such an equation 
the leading order coefficient, and generically only this one, transforms as a tensor, if the entire 
equation does [s^. This will render definition ([2]) below covariant. For the example of the geometry 
being given by an inverse metric tensor field (G"^ = G^'^ with non-zero determinant everywhere) 
and a scalar field (j) {d\mV = 1) obeying the massless Klein Gordon equation G'^^"''^da^da24' " 
^"'^mnG'^'^da^(j) = 0, where T are the Christoffel symbols of the metric Gab, we indeed find that the 
leading quadratic order coefficient transforms as a tensor, while the linear order coefficient does 
not (and could not, since it must ensure that the enrire equation transforms as a scalar). 

Requiring that matter equations of the form ([T| are predictive, interpretable and quantizable 
imposes necessary conditions on the underlying geometry G. These conditions have been derived 
and explained in detail in Q]. Here we present a practical summary of these conditions and their 
implications as far as they are directly relevant for the present article. All constructions revolve 
around the totally symmetric contravariant tensor field P defined from the leading order coefficients 
of the matter field equations ([I]) by 

P'^ -'^^^^{x)h, . . . h,^^, := pdet {Q[G]%-^^{x)ka, . . . K^) (2) 

for all points x £ M and cotangent vectors k G T*M and a scalar density function p constructed 
from the geometry G such as to be of opposite density weight to the determinant over the tensor 
representation indices. To lighten the notation, we will often use the shorthand P{x, k) for the left 
hand side of Eq. ([2]). Furthermore we may agree, since no information is lost and it is technically 
convenient, to remove any repeated factors into which the field P may factorize; so if the above 
construction yields P{x, k) = Pi{x, k)°'^ ■ ■ ■ Pf{x, k)"f then we consider instead the reduced tensor 
field P defined by P{x, k) = Pi{x, k) ■ ■ ■ Pf{x, k). The physical meaning of the tensor field P is 
revealed by the eikonal equation [20] for the dynamics ([T]), which shows that 

P{x,k) = (3) 

is the dispersion relation that a covector k £ T*M must satisfy in order to qualify as a massless 
momentum. For our simple example of a Klein-Gordon field on a metric geometry, the determinant 
in ([2]) is of weight zero, and for the choice p = 1 we obtain P*i*2 = G*^*^, and one indeed recovers 
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the familiar massless dispersion relation C^^"^^ kaj^ka2 = 0. An instructive non-metric example is 
provided by abelian gauge theory coupled to an inverse area metric tensor geometry [ll|, [3] , which 
is based on a fourth rank contravariant tensor field G featuring the algebraic symmetries G"'^'^'^ = 
Qcdab _ _Qbacd.^ calculation of the principal polynomial (after removing gauge-invariance, obser ving 



Id, 



resulting constraints on initial conditions and re-covariantizing the expression) one obtains 
\\\ in dim M = d dimensions the totally symmetric tensor field 

ph---i2(d-2) = p(^Q'j ^^^^ ^g^^ ^ ^Qaa.ibi{iiQi2\a2b2\i3 ^ ^ ^ Qi2(d-3)\a-d-2bd-2\i2d-5Qi2(d^2))<^d-lbd-lb 

(4) 

of tensor rank degP = 2(d — 2), with some scalar density p{G) of weight +2 constructed from G. 
In four spacetime dimensions, for example, where the area metric may be decomposed into a cyclic 
part Gc with Q'^^'^^ = q and a totally antisymmetric part given in terms of a scalar density ^ of 
weight -1, G"^^'^ = G^*^'^ + ^e"^""^, one may chose p{G) = -l/(24$2). This non-trivial example 
for a field P illustrates two salient points. First, it reveals what a dramatic accident it is that in 
Lorentzian geometry the field P, which will be central to all further developments, is essentially 
identical to the fundamental geometric field g; for in area metric geometry, one not only sees that 
P is a tensor field of generically entirely different tensor rank than the underlying fundamental 
geometric tensor G, but may also feature an entirely different index symmetry structure: the tensor 
P is always totally symmetric. Second, it exemplifies the rule that generically the fundamental 
geometry G cannot be reconstructed from the field P {sil]. With these remarks on the role of the 
field P as a dispersion relation and its origin in matter field equations coupled to some tensorial 
geometry, we are now prepared to lay down the three crucial algebraic conditions that the tensor 
ield P needs to satisfy. These conditions in turn restrict the geometric tensor G that underlies P 



The first condition, predictivity of the matter field equations, translates into the algebraic 
requirement that the tensor field P be hyperbolic {isl. IigI. This means that there exists a covector 
field h with P{h) > such that for all covector fields r there are only real functions A on M such 
that 

P{x,rix) + X{x)h{x))=0 (5) 

everywhere. Obviously if P{h) < 0, one could arrange for P{h) > simply by changing the overall 
sign of the density p appearing in ([2]), and we will agree to do so for definiteness '33!]. In any 
case, it is useful terminology to call a covector field h, if it indeed exists, a hyperbolic covector 
field with respect to P. Only hypersurfaces whose canonical normal covector fields (defined to 
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annihilate any tangent vector field to the hypersurface) are hyperbolic can serve as initial data 
surfaces for equations of the type ([I]). We will return to this point later. For our example of a 
metric geometry, it is easy to check that P*i*2 — (7«i«2 hyperbolic if and only if the inverse metric 

has Lorentzian signature (H ■ ■ — )) and that the hyperbolic covectors are exactly those covectors 

for which P^^^^hi-^hi^ > 0; in other words, initial data surfaces need to be spacelike. The reader 
be warned, however, that such a simple characterization of hyperbolic covectors and thus initial 
data surfaces merely by the sign of their co-normals or tangent vectors under P is not generic and 
merely a coincidence in the metric case. The underlying general definitions however work for all 
geometries. 

The second condition, interpretability of the matter field equations translates into a time- 
orientability condition for the underlying geometry. This is simply the algebraic requirement that 
also the so-called dual tensor field be hyperbolic. Indeed, for any hyperbolic tensor field P, 
one can show that there always exists a totally symmetric covariant dual tensor field of some 
rank degP* such that 

P*{x,DP{x,k{x))) =0 (6) 

for all covector fields k with P(x,k{x)) = everywhere, where DP{x,q) denotes the vector with 
components DP{x, /c)" = {degP)P'^^'^'"^'^°^Pkf^ . . . kj^ and we used a shorthand for the evaluation 
of the field on a vector that is analogous to the previous one for P on a covector. The dual 
tensor field P"^ is unique up to a real conformal factor, and can always be constructed, essentially 
by determining a Grobner basis {l?!. For our example of metric geometry, a dual of P*2*2 = (7*1*2 
is given by Pj[i,^ = Gi-^i^, as one easily verifies. Returning to the general case, time-orientability 
means that there exists a vector field H such that for every vector field R there are only real 
functions fi on M such that 

P*{x,R{x) + n{x)H{x)) = (7) 

everywhere. A vector field H satisfying this condition will be called a time-orientation. Once a 
time-orientation has been chosen, it is useful to consider, separately in each tangent space, the 
connected set Cx of all hyperbolic vectors to which the vector Hx of the time-orientation belongs. 
According to a classic theorem [l^, constitutes an open and convex cone in the tangent space 
TxM, and we will call Cx the cone of observer tangents (to observer worldlines through the point 
x). Note that in general, hyperbolicity of P does not already imply hyperbolicity of P"^, and 
thus predictivity does not imply time-orientability in general. For our metric example, however, 
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it trivially does; Pj*^ = Gi-^i^ is hyperbolic if and only if P*i*2 = G^^^^ ig^ and the cones C* of 
observer tangents are the timelike vectors X at each point that are future-oriented with respect to 
some global timelike vector field T representing the time-orientation, i.e., satisfy Ga^aj^"^^"^ > Oj 
again, this simple sign condition to decide membership of X in the cone of observer tangents 
selected by T is a coincidence in the metric case, and again has to be replaced by the underlying 
general definition above for other geometries. 

The third condition on the matter field equations, namely that these are quantizable, is that the 
geometry be energy-distinguishing. This simply means that all observers agree on the sign of the 
energy of a massless momentum. More precisely, a geometry is energy-distinguishing if for every 
point X G M and every massless momentum k either k{X) > or —k{X) > for all X G Cx- In 
a hyperbolic, time-orientable and energy-distinguishing geometry, one can then also meaningfully 
define massive particle momenta of positive energy at some point x as those hyperbolic covectors 
q G T*M for which q{X) > for all observers X G Cx- To be able to do this is of crucial 
importance when performing a split of a basis of solutions to the field equations into positive and 
negative frequency solutions in a canonical quantization of the matter field. These massive positive 
energy momenta constitute an open convex cone Cx in the cotangent space at x. The mass m of 
such a positive energy massive particle momentum q £ Cx is then defined by 

P{x,q) =m'^^^^ . (8) 

In Lorentzian metric geometry, the above definitions of course recover as the positive energy massive 
and massless momenta precisely those timelike and null covectors whose application to a future- 
directed timelike vector is positive. It may be worth emphasizing again that for a covector to qualify 
as massive, it must not only satisfy the massive dispersion relation ([8]) but indeed be hyperbolic, as 
stipulated above. Only in Lorentzian geometry does the massive dispersion relation already imply 
hyperbolicity. 

Only if a geometry satisfies the three conditions laid out above can one associate worldlines 
with the massless and massive dispersion relations. For only then can one solve for the momenta 
q after variation of the Hamiltonian actions 

Snia.ss\ess[x,q,X\ = J dr [x°-qa - XPiq)] and Sniassivc[x,q, X] = j dr x°-qa - Xln P{-^) , 

(9) 

respectively. Defining the Legendre map Lx for all covectors q in the open convex cone Cx of 
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positive energy massive momenta at some point x as 

L'^ix, q) = ^(^'^'^ , (10) 

which by virtue of the energy-orientabihty of P possesses a unique inverse on its domain, 
one finds jfl] that the worldlines of free massless and massive particles are stationary curves of the 
reparametrization-invariant Lagrangian actions 

Smass\ess[f^,x] = j dr nP*{x) and Smassivc[x] = m j dr P{L'^{x)y^^ , (11) 

respectively. The massive particle action reveals the physical meaning of the Legendre map L, since 
one readily derives that the canonical momentum of a positive energy massive particle is related 
to the worldline tangent vector as q = mL^^{x) if one chooses the proper time parametrisation 
P(L-i(x)) = 1 along the worldline. Put simply, the Legendre map raises the index on a positive 
energy massive momentum, in one-to-one but non- linear fashion. For the example of Lorentzian 
geometry, we find that under the familiar proper time parametrisation Gabx"'x^ = 1) the worldline 
tangent vector x and the corresponding particle momentum q of a particle of mass m are related 
through mx"' = C^^qi,. The massless and massive Lagrangian actions for the free point particle on 
a Lorentzian spacetime recover the standard textbook postulates. 

Of central importance for the aim of this article, namely to derive the dynamics of hyperbolic, 
time-orientable and energy-orientable geometries, is the following insight. Hypersurfaces that are 
potential carriers of initial data and at the same time accessible by observers are those whose co- 
normal at each point lie in the cone L~^{C'^). This is because, on the one hand, the purely spatial 
directions seen by an observer with tangent vector X G are precisely those vectors annihilated 
by the covector L^^[X). On the other hand, the cone of these observer co-tangents can be shown 
to always lie within the cone C of hyperbolic covectors for hyperbolic, time-orientable and energy- 
distinguishing geometries. In Lorentzian geometry, such inital data surfaces accessible to observers 
are simply the spacelike hypersurfaces. Incidentally, only when L~^(C*) does not only lie within 
C, but entirely coincides with it, is the theory free of particles travelling faster than the speed of 
some light 34i |. 



The deformation of such observer-accessible initial data hypersurfaces, separately in normal and 
tangential directions, will be the topic of the next section. While a generic hypersurface directly 
gives rise to tangent directions, but merely normal co-directions, it is only the Legendre map (and 
thus the spacetime geometry) that allows to associate a normal co-direction n of a hypersurface 
with a corresponding normal direction T = L{n) if n lies in L^^(C'^). Normalizing the latter 
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by requiring P{L~^{T)) = 1 corresponds to requiring that the normal direction be tangent to an 
observer worldhne with proper time parametrization. Thus the normal deformation of observer- 
accessible initial data hypersurfaces feels the spacetime geometry through the Legendre map. It is 
this role of the Legendre map that we will see to hold the key to the derivation of the gravitational 
dynamics for general tensorial spacetimes. 



B. Deformation of initial data surfaces 



The aim of this paper is to find dynamics that develop initial geometric data from one initial data 
hypersurface to another, such that sweeping out the spacetime manifold in this way one reconstructs 
a hyperbolic, time-orientable and energy-distinguishing dispersion relation everywhere. In this 
section we describe initial data hypersurfaces by their embedding maps and study how functionals 
of this embedding map change under normal and tangential deformations of the hypersurface. The 
functionals of interest later on will be the induced geometry seen by point particles in section IIIII 
or the induced geometry seen by fields in section IIVI The change of a generic functional of the 
embedding map can be expressed by a linear action of deformation operators on such functionals, 
and it is the commutation algebra of these deformation operators that we are after js^. 

More precisely, we consider a hypersurface X{T,) defined by an embedding map X : T, ^ M ol a 
smooth manifold S with local coordinates {y"} into the smooth manifold M with local coordinates 
{x'^}; here and in the remainder of this paper, latin 'spacetime' indices run from to dimM — 1 
while greek 'hypersurface' indices run from 1 to dimM — 1. Without additional structure, the 
embedding defines at each point y of the hypersurface dimM — 1 spacetime vectors 

tangent to the hypersurface X^E), which in turn define, up to scale, normal spacetime covectors 
n(y) as the annihilators of all tangent vectors, 

n{y){ea{y)) = a = 1, . . . , dimM-1 . (13) 

Only if we restrict attention to initial data hypersurfaces whose data are accessible to observers, 
by requiring that the n{y) lie in the respective cones -L~^(C*) everywhere along the hypersurface 
can we impose the normalization P{n{y)) = 1 and thus obtain a unique spactime vector 
field T{y) = L{n{y)) representing the normal directions, rather than normal co-directions, away 
from the hypersurface. Thus an accessible initial data hypersurface X{Ti) induces a complete 
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spacetime tangent space basis {T{y), ei(y), . . . , edimM-i(y)} at every of its points, and dual basis 
{n(y), e^(y), . . . , e'^^"^ '^^^^ (y)} in cotangent space. 

We now consider deformations of the hypersurface X{T,). Technically, this is done by prescribing 
a smooth one-parameter family Xt of embedding maps such that the original embedding map 
X is recovered for t = 0. Then the connecting vector field dXf/dt in spacetime, between the 
hypersurfaces XfiY^) of this family, can be uniquely decomposed along the undeformed hypersurface 
into a sum of a purely spatial and a purely tangential part, 

X(y) = iV(y)r"(y) + iV°(y)eS(y), (14) 

where the hypersurface scalar field N and hypersurface vector field components N'^ are given by 

N{y)=n{y){X{y)) and iV"(y) = e°(y)(X) (15) 

and thus completely parametrize any small deformation of the embedding map X into X + dtX . 

The linear change of functionals under changes of the embedding map is conveniently studied 
in terms of normal and tangential deformation operators. More precisely, we define the normal 
deformation operator 

n{N) = JjyN{y)T%y)jj^^ , (16) 

acting on arbitrary functionals F of the embedding function. The change of such F under the 
deformation ()14p is then given to first order hj'H{N)F. Similarly one obtains for a purely tangential 
deformation the first order change 'D{N°'da)F through the tangential deformation operator 

P(iV-9„) = jjyN%y)el{y)j^^ . (17) 

A trivial check on the geometric meaning, which this terminology attaches to these operators, 
is their action on the components of the embedding map itself; with the definitions of the delta 
distribution and functional derivatives one finds 

'H{N)X''{z) = N{z)T''{z) and P(iV°9c.)^"(^) = iV°(2)eS(^) , (18) 

which indeed are precisely the normal and tangential components of the deformation (|14p . Since 
the embedding is a linear functional of itself, this shows that ()16p and ()17p indeed are the operators 
that bring about the normal and tangential deformations of functionals to linear order, as desired. 

Finally we may calculate their commutator algebra. The latter will play a crucial role througout 
this paper. Now the basis vectors T and ei, . . . , edimM-i are functionals of the embedding map, 
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and it is thus clear that multiple application of deformation operators will require to know their 
functional derivatives with respect to the embedding functions. While for the hypersurface tangent 
vectors one obtains 



8Xb{z) 



(19) 



in straightforward fashion directly from their definition, one needs to work somewhat harder from 
the definition of T to find 



6Xb{z) 



(degP-l)(e>,P"'')(y)9;35,(z) 



+ 



degP-1 
deg P 



{y)6y{z), (20) 



where in the first summand one of the hypersurface tensors defined in (|34p appears. Note that the 
dispersion relation enters only in the variation of the normal vector, but not of the tangent vectors. 
This is because the definition of the former employs the Legendre map defined by the dispersion 
relation. This is indeed the way t he g eometry enters into the deformation algebra, which is now 
straightforwardly calculated to be [36 1 



[n{N),niM)] = -P((degP-l)P"^(M9/3iV - NdpM)da) , (21) 
[V{N''da),'H{M)] = -niN^'daM) , (22) 
[D{N'^da),V{M'^dl3)] = -V{{Nl^di3M'^ - M^dpm)da) ■ (23) 

The exclusive appearance of the hypersurface tensor field components P°'^ induced from the space- 
time tensor P by virtue of P°'^ = P'^^^^ '^'^^sP-^e^e^nj-^ . . ■nf^^^p_^ on the right hand side of the 
commutator of two normal deformation operators — the last two commutators are indeed fully in- 
dependent of the hypersurface geometry — originates entirely in the use of the Legendre map when 
defining the spatial fields. Thus P"^ appears irrespective of which type of geometry on the hy- 
persurfaces one chooses to study (possible choices are the pullbacks of the geometry seen by point 
particles considered in chapter IIIII or the geometry seen by fields considered in chapter lIVp . The 
calculation of the on the right hand side of the algebra above in terms of the hypersurface 
geometry is just more complicated for geometries seen by fields than for geometries seen only by 
point particles, but it is always the P"^^ that appears there. Finally it is useful to observe that the 
tangential deformation operators constitute a subalgebra. 
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C. Towards canonical dynamics for hypersurface geometries 

So far in this paper, we tacitly assumed to have knowledge about the values of the geometric 
tensor G — and thus also the cotangent bundle function P derived from it — at every point of the 
entire spacetime manifold M. This enabled us to derive how any functional F of a hypersurface 
embedding map X : ^ M changes under a change (jl4p of the embedding map. We are 
particularly interested in the particular type of functionals of the embedding map that arise as 
normal and tangential projections of a spacetime (r, s)-tensor field F to an embedded hypersurface. 
For simplicity, consider a (l,0)-tensor field F on M, which induces the projections 

F^{y)[X]:= F{n{y)) and F°(y)[X] := F(e"(y)) for a = 1, . . . , dimM - 1 , (24) 

which yields the collection of functionals F^{y) = {F'^ (y) , F'^ (y)) , where we used the spacetime 
covector frame {n, e^, . . . gdimAf-ij g^jQj^g i^j^g hypersurface to project F. One proceeds analogously 
for tensor fields of valence {r,s). Knowing the value of the tensor field F throughout spacetime, 
and in particular in a neighborhood of an embedded hypersurface X{T,), we can write the linear 
change of the functionals F^(y) under a deformation of the original hypersurface controlled by the 
lapse N and shift A^" as 

/ dz [N{z)n{z) + iV"(z)P,(z)] F^{y)[X] , (25) 

where we introduced the localized operators T-l{z) := Tiidz) and Va{z) := V{6zda)- 

But this omniscient view of the values of the tensor field F and the geometry G at every point of 
spacetime is not afforded by us mere mortals. What we have access to, at best, are the values F^{y) 
on S, understood as mere hypersurface tensor fields, rather than functionals of the embedding map. 
If we then wish to predict the values of the F"^ on some different hypersurface through spacetime, 
we need to stipulate how these fields change from the initial hypersurface X{Ti) to a deformed one 
near-by, and we will see in a moment that we are rather constrained in the way we can stipulate 
such equations of motion. Anyway, since we are ignorant of any of the field values of F away 
from the hypersurface, we need to compensate for this lack of knowledge by adjoining canonical 
momentum densities (pA of weight one to each configuration variable -F"^, which is equivalent to 
introducing a Poisson bracket 

{C,b] := [ dz 
is 

on the space of functionals of the phase space variables {F^,(f>A), which is sometimes referred to 
as superspace. One can then give dynamics to the hypersurface fields F^ by stipulating that their 



6C 6D SD 6C 



6F^ 



6F^ 5(1) A 



(26) 
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values change by the amount 

{F^{y), I dz N{z)'H{z) + N"{z)V^{z) ] (27) 



s 

when evolved to a neighboring hypersurface whose deformation from the initial one is determined 
by the lapse N and shift A^", where the quantities % and are some a priori arbitrary func- 
tionals of the phase space variables [F'^,cj)A)- For brevity, and in accordance with the standard 
terminology in geometrodynamics, we will refer to T-L{y) as the superhamiltonian and to T>a{y) as 
the supermomentum. The dynamics ()27p are further assumed to be supplemented by first class 
constraints 

?i(y)«0 and P„(y) « (28) 

implementing the required diffeomorphism gauge symmetry. 

It is clear that if the dynamically evolved hypersurface field values on the deformed hypersurface 
are to coincide with what the hypersurface deformation (I25|) yields, independent of any particular 
deformation {N,N"'), then we must require that (j25p coincides with (127p . or equivalently, 

n{y)F%)[X] = {F%),niz)iy)} and V^iy)F%)[X] = {F%),V^{z){y)} . (29) 

We cannot extend these equations to the momentum variables, since we do not know at this stage 
how the latter can be understood as functionals of the embedding map — this is for the dynamics to 
determine. But using the relations ()29p in the deformation algebra ()2ip . p2p . ([23]) and the Jacobi 
identity for the Poisson bracket, one finds that a sufficient condition for our above compatibility 
requirement is that the functionals 

n{N):= f dyN{y)n{y) and P(iV"9„) := / dyiV°(y)P„(y) (30) 

of the phase space variables have Poisson brackets that represent the deformation algebra commu- 
tation relations [s^] 

{n{N),n{M)} = V{{degP-l)P''P{MdpN - NdpM)da) , (31) 
{V{N'^da),n{M)} = niN'^daM) , (32) 
{'D{N'^da),V{M^dp)} = ViiN^dpW - M^dpN'')do) . (33) 

The extent to which this representation requirement is not necessary to satisfy our compatibility 
condition, however, precisely encodes the information concerning the functional dependence of Ti 
and V on F^, while ([29]) already determines 5'H{z) / ScpAiy) and S'Da{z)/S4>A{y)- We will return to 
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this point when constructing % and T> from the above algebra, and indeed the major part of the 
remainder of this paper wih be devoted to this construction. 

At this point the paper sphts into two different projects. While both are concerned with finding 
gravitational dynamics — by way of finding the supermomentum and % satisfying the Poisson 
algebra (f3T]l . (p2]) . ([33]l according to the programme laid out in this section — they do so for different 
geometric degrees of freedom. The first project, treated in chapter IlIU considers the projections 
of the components of the dispersion relation P as the geometric degrees of freedom on a hyper- 
surface, and derives the equations that determine its superhamiltonian and supermomentum. In 
contrast, the second project, treated in chapter \IV\ illustrates how the projections of the funda- 
mental geometric tensor G are taken as the degrees of freedom to which one gives dynamics. This 
point of view is more fundamental, but its details depend heavily on the algebraic properties of 
the geometric tensor. 
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III. DYNAMICS OF DISPERSION RELATIONS 

Point particles only see those aspects of a tensorial spacetime geometry G that are encoded in the 
tensor field P underlying the massive and massless dispersion relations. It is a mere coincidence in 
Lorentzian geometry that the tensor field — Qaia2 contains precisely the same information as 

the fundamental geometric tensor field Gab to which fields couple. For any other tensorial spacetime 
geometry, the totally symmetric tensor field P can be expressed in terms of the fundamental 
geometric tensor field G, but not vice versa. So if one is interested in a full gravitational theory 
to which fields and point particles can couple, one needs to derive dynamics for the fundamental 
geometric tensor G, and we will do so in chapter IIVI But if one is only interested in gravitational 
fields acting on, and generated by, point particles, one may alternatively construct gravitational 
dynamics directly for any hyperbolic, time-orientable and energy-distinguishing tensor field P. 
Indeed, we obtain a rather sweeping result in this chapter: we derive a system of homogenous 
linear partial differential equations whose solutions yield all possible canonical dynamics for physical 
dispersion relations. 

A. Phase space for geometries seen by point particles 

Since we wish to study the spatial geometry seen by point particles on an initial data hypersur- 
face X{T,), we are interested in the hypersurface tensor fields P, P'^^ , . . . , p^i^a- .Odcgp that arise 
as functionals of the embedding map through 

pa,...a, p^^a, ^ ^ ^ ^ ^ ^a, ^(y)^ ^(y)) for / = 0, . . . , dcg P , (34) 

^ V ' ^ V ' 

/ dcg P-I 

where we used the complete spacetime covector frame {n, e^, . . . , gdeg Af-i j g^jQj^g the hypersurface 
to project the spacetime tensor field P'^i---'^dogp onto S. The normalisation conditions P{n) = 1 
and T"e" = immediately imply that the two functionals P and are constant along S. This 
property is conserved under hypersurface deformations and thus the P and P°^ do not carry any 
dynamical information and can be discarded as configuration variables. The remaining hypersurface 
tensor fields however allow, in combination with the dual basis, the complete reconstruction of the 
spacetime dispersion relation at every point of the hypersurface and thus, in their entirety, present 
the geometry on the hypersurface — as it is seen by point particles. 

As we have laid out in section III C\ the point of dynamics is to generate, not assume, the 
values of P throughout the spacetime manifold, starting only from initial data in form of totally 
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symmetric tensor fields 

paia2 pa\a2a-i paia2---adogP (35) 

on a manifold T, of dimension dimM — 1. To lighten the notation, we collect the fields (j35p 



in a quantity carrying a multi-index A = (aia2; aia2a3; aia2 • • • Odcgp) consisting of 
totally symmetric tuples of spacetime indices. Starting from these configuration variables, which 
correspond to initial data on only one hypersurface Xo{T,) say, the lack of knowledge about the 
value of the configuration variables on neighboring hypersurfaces Xt{T,) must be compensated for 
by adjoining canonical momenta 

or TTA for short, to the configuration variables (j35p on S. This introduction of canonical momenta 
is of course equivalent to adopting a Poisson bracket 

6F 5G 6G 5F 



{F,G}= / dy 



(37) 



on the space of functionals of P"^ and tta, whose elements we denote with a hat. The configuration 
variables P"^ and momenta tta are trivally functionals of themselves, and this is the reason why we 
made them carry a hat from the beginning. In order for the integral (j37p to be well-defined under 
changes of chart on S, the momenta must be tensor densities of weight one. This ultimately follows 
from the definition of the functional derivative. We now set out to determine the supermomentum 
V and the superhamiltonian H that close according to (fSTI) . (f32ll and (f33ll and evolve the canonical 
variables (P"^,7r^). 

B. Construction of the supermomentum 

The strategy to determine the superhamiltonian and supermomentum from the Poisson 
algebra — in accordance with the philosophy laid out in section III CI — follows from the fact that 
the supermomentum functionals constitute a subalgebra that must be solved without recourse to 
the further relations. This is what we do in this section. More precisely, using the definition ()17p 
of the tangential deformation operator we may calculate the change of the functionals (j34p under 
tangential deformations. In order to do that, we need to know the functional derivatives of the 
dual hypersurface basis vectors {n, e"} with respect to the embbeding map. For the hypersurface 
conormal, we have 

Sna{y) _ 1 



6Xb{z) degP 



{nan,, . . . n,-, ^a6P^i-^<^-^)(y)5,(z) + miy)e2{y)djy{z), (38) 
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while the functional derivatives of the e" read 

- {naeln,, . . . n,^^^^d,P^--^'^'^^n[y)5y{z) . (39) 

Using these results on the left hand side in the second equation of (p9|) . and the second weak 
equality in ([28l) . one calculates that 

{p%),V{N'^d^)} = {C^P)^{y). (40) 

This is entirely plausible, since the deformation operator, and thus its represention in form of the 
supermomentum, push the configuration variable along the hypersurface vector field N, and this 
is precisely what the Lie derivative is defined to do. But then it follows from the Jacobi identity 
for the Poisson bracket ()37p . followed by a functional integration with respect to vr, that 

{TTAiy), ViN'^d^)] = {Cj^7f)A{y) , (41) 

for the covariant tensor densities tt^. Again this is more than plausible, since the tangential 
deformation operator merely reshuffles all the initial data. In summary, we obtain a set of pairwise 
coupled functional differential equations for all / = 2, . . . , deg P, namely 
6V{N) 



^pQl...Q/(y) 

6V{N) 



{-d/,N^7:a,...aj - N^d)37ra,...aj - Id^^,NPK,...o.My) , (42) 

pai...«i _ J5 Ar{"ip°2...a/)/3)(y)^ (43) 
OT^ai...ai{y) 

which are integrable since all second functional derivatives of T){N) commute. These equations are 
uniquely solved by 

degP 

i>{^)= / (iyA^^(y) [5;3P"i - "^7r,,...,, +I5«,(P"i----^7r„,...,,^)] , (44) 
1=2 -'^ 

where an a priori non-zero additive integration constant is forced to be zero by (jSSp . This is already 
the desired supermomentum appearing in the dynamics (I27p . Note that, in the case of degP = 2, 
the supermomentum reduces to the standard form = 2P^'^'S/ isTta'y known from general relativity 



C. Construction of the non-local superhamiltonian 

Now that the supermomentum, and thus the right hand side of the bracket ()3ip is known, we 
can start to determine the superhamiltonian by extracting the information contained in this and 



21 



the other remaining bracket (j32p . The latter simply tells us that the superhamiltonian must be a 
scalar density of weight one. For on the one hand, we concluded from (I40p and (I4ip that {■,'D{N)} 
acts like a Lie derivative on any functional of the geometric phase space variables, and thus on the 
superhamiltonian in particular. But on the other hand, letting B{z) = 6y{z) in the bracket (j32p 
we obtain 

{n{y)MN)} = d^{niy)N^{y)) (45) 

after an integration by parts on the right hand side. But this is the Lie derivative of H only if the 
latter is a scalar density of weight one, as claimed, and this is all that follows from this second 
bracket. 

Again we approach the solution of the Poisson relations starting from one of the necessary 
relations (|29p . From the left equation we indeed find that 

(/-degP)P^"i-"^(z) + (degP-l)/p("2...a7(^)pai)/3(^)l ^ 

(46) 

where the expression in [. . . ] contains the configuration variables P^, all frame covectors, and their 
first derivatives. The second term, in contrast, only contains the configuration variables, and this 
will become important shortly. If the second bracket vanished in general (which, however, only 
is the case for degP = 2), the superhamiltonian would be a function, rather than a functional, 
in the momenta tta according to (I46p . However, the simple form of the djsN term allows to 
directly determine the non-local contribution of the momenta to the superhamiltonian. In fact, it 
is straightforward to check that one can decompose the superhamiltonian as 

n{y)[P, vr] = ?ilocal(y)[^](^) + -Wnon-local(y)(^, Ott) , (47) 

namely into a local part ?^iocai(y)) which is indeed a functional of P but only a function of tt, and 
the explicit non-local part 

deg P 

^non-iocai(y)[^',vr] = ^ [(deg P-/)9^ (P^"^ •••°^7r„, ,,,,,)- (deg P-1) / 9^ (P^^-"^ P"i^7r«, ...«,)] (y) , 

1=2 

(48) 

which is thus a completely known functional of P and tt that generates the non-local second term 
in (j46p . Note that the non-local part 'Hnon-local(y) of the superhamiltonian is the divergence of a 
vector density of weight one and thus a scalar density of the same weight. Hence, the decomposition 
(j47p turns the superhamiltonian into the sum of two tensor densities of weight one [l9]. This means 
that we reduced the problem of finding the superhamiltonian as a functional of both phase space 



OH(I\ ) 



Niz) 



+ dpN{z) 



22 



variables P and vr to the much simpler problem, as it will turn out, of determining the local part 
that is a functional in P but only a function in tt. In particular, this will allow to make a power 
series ansatz for T^iocai iii ^r. 



D. Lagrangian reformulation 

At this point we explicitly know the supermomentum D and the non-local part T^non-iocai of 
the superhamiltonian (j47p . The still undetermined local part 7^ local of the latter enters the only 
remaining Poisson bracket (j3ip quadratically on its left hand side, 



dz 



s 



'^'^(a;)local _^ (5?^(x)non-local 

5P^{z) SPa{z) 



5'H{y)\oca\ ^ non-local 



(x o y) , (49) 



5tta{z) 5t^a{z) 

where the contributions of the non-local part of the Hamiltonian are explicitly known from taking 
the functional derivative of (j48p . Here and in the remainder of this paper, repeated multi-indices 
indicate sums of the form 

degP 

C'^Da = • (50) 

1=2 



The quadratic appearance of ^locai in ([Ml) seriously complicates a power series ansatz for it in the 
momenta vr. Remarkably, a Legendre transformation [4^ replacing the momenta t:a by Legendre 
dual variables 

d-KA^x) 

from which conversely the momenta depend as a function, 7r(x)[P](X), allows to turn the equation 
l9j) that is quadratic in T^iocai into an equation that is linear in the "Lagrangian" 



L{x)[P]{K) := 7r^(x)[P](K) K^{x) - n{x\,UP] {^[P]{K)) , (52) 
since then one finds 



local 



^"^^"^ J#, = vr^(x)[P](K). (53) 



5PA{y) 

Let us further define the coefficients Qa^^ and M^'' by 

^^non-local(x) _ (^^^^^^ (^^^^ _ , (54) 



SP'^iz) 

= M^«(y)9c<5,(.) - d^M^^{y)5y{z) , (55) 
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which yields the expressions 



^; Qi ...i 



5f+i(degP-/)5ff--^^ 



, -5fl(degP-l)p('^2...A^/^iV 
-6f_,I{degP-l)P^^^^5^,l::i'l^ , (56) 



depending only on the configuration variables P. Rewriting the Poisson bracket (|49p with the help 
of the Lagrangian and integrating out the appearing delta distributions, its left hand side takes 
the form 



5L{x) 



5L{x) 



+ M^Hy)QA''^{x)d^7rB{x)d^6y{x) 



- K^{y)QA''^{x)TiB{x)dpS^{y) + Qa""^ {x)T^B{x)M^^{y)dl^5y{x) 
+ QA''^{x)7:B{x)d^M''^{y)dp5^{y) - (x o y) , 



while the right hand side becomes (degP — 1) times 



degP 



1=2 



{y)daSx{y) - (x^ 



where vr is given by the second of the equations (j53p . A key observation is now that the dependence 
of the terms in square brackets on the right hand side may be changed from y to x while the 
dependence of the delta distribution multiplying it remains unchanged; due to the exchange term 
(x o y), the resulting distributions are the same. The same remark applies to changing the 
dependence of Qa^^ and 971"^ from x to y in the third term on the left hand side. We may thus 
collect the derivative terms dn from both sides into an expression of the form 



T^^'''ix)iP)df,6,{y)d,TrAix) - (x o y) 



(58) 



on the left hand side of the original Poisson bracket relation. Crucially, one finds that T'^'"' = T^"'' 
by inspecting the explicit expression for the above coefficients. It is only due to this fact that ([58|) 
is equal to 



T^^''{x){P)d^dJ,iy) - d^T^^-'ixmd^d^iy)] 7:a{x) - (x o y) 



(59) 



as a distribution in two variables. Thus all dn terms can be made into local expressions in the K 
by virtue of the second relation in ([53]) . 
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Combining all terms of the original Poisson bracket (j3ip in this fashion, one obtains its entirely 
equivalent formulation as a homogeneous linear functional differential equation in L = L[P]{K), 



= -^i^K^{y)+d^c 



where the coefficients [/^'"^ contain the configuration variables, 

whereas the coefficients S^'^ also contain their first partial derivatives, 

gm...ar,. = -(degP- f)P''^a^P"i-"^ +/(degP-/)(degP- f)P^("i-a^P"^)'' 

+2/(degP - + /(degP - l)(degP - 2)P'^''("ia;3P"2...a/) 

-/(/ - l)(degP - l)2p/^("ip"2...a/-i^^pa/)/3^ (^Q2) 

where in the case degP = 2 the last term is to be read as —2P^^"^dpP'^'^^^. Once one has solved 
()60p for the Lagrangian L[P][K), one can recover the momenta 

„ dL{x)[P]{K) 

^a{x) = g^A(^) ^ (63) 

conversely expressing K^{x) = K^{x)[P]{7r), and one also finds the local part of the superhamil- 
tonian as 

^(x)iocai[^](vr) = nA{x)K^{x)[Pm - L(x)[P](i^^(x)[P](vr)) . (64) 

This then of course amounts to the full determination of the gravitational dynamics, since the 
supermomentum and non-local part of the superhamiltonian are already known from previous sec- 
tions. But the difficulty of solving (I60p consists in this being a distributional functional differential 
equation for L. 



E. Reduction to differential equations 

In this section we will reduce the equation (j6Up to a countable set of linear partial differential 
equations for the functional L that determines the still missing local part of the superhamiltonian. 
This reduction takes place in two steps: 
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distributional functional differential equation 
i 

distributional differential equations 

i 

differential equations 

The first step to achieve this exploits the linear homogeneous structure of the equation (|60p by 
making a power series ansatz 

oo 

L{x)[P]{K) = ^C{x)[PW..A^K^Hx)...K^^{x) (65) 

i=0 

with coefficients that are so far undetermined functionals of, and this is the essential point, only 
the configuration variables. A power series expansion is justified since we took care in constructing 
L as a mere function of K, while it remains a functional of P. Since the velocities are defined 
as the partial derivatives of the weight-one scalar density Hiocai with respect to the tensor densities 

of the same weight, the velocities themselves are tensors. The coefficient functionals C[P]ai...An 
are thus tensor densities of weight one just as the Lagrangian L[P]{K). 

Insertion of (|65]l into (f60]l replaces the latter, a distributional differential equation for L[P]{K), 
by a countable set of such equations for the coefficient functionals C[P]ai...Ajs,'^ one equation for 
each order in i^. Extracting the A-th order equation by application of the functional derivative 
operator 

5KBi{xi)...6KBm(^xn) ^^^^ 

to the equation (I60p and evaluating the result at A' = 0, we will now see that one obtains at 
A-th order a distributional equation in A + 2 variables x,y,xi, . . . ,X]sr. Indeed, the zeroth order 
contribution is 

6C{x) 



= dyc 



6P^{y) 

while the contribution at order A > 1 is 



+ C{x)a [U^^''{x)d^d,.6,{y) + - (x o y) , (67) 



= {(AT + 1)IG{x)ab,...b^ (f/^^"(x)a^5,5,(y) + S^^'{x)d^5,{y)) 



+ A! dyc 



\^AC^y^^C{x)B,...B, 



N 6C{x)^^ 



5P^{y) 

(y) 



NmQ^B,^'^ix)CXx)B,...B^)MdpS,{y)} 
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where ~ instructs to omit a term. As usual, these distributional equations are to be understood 
by first applying them to test functions f{x, y,xi, . . . ,xn) and then integrating over all variables. 
In order to convert the thus constructed functional differential equations into regular differential 
equations, we restrict attention to coefficients C{x)ai...aA^] that are determined by the value of 
P and all its derivatives at x, so that 

c{x)[PW..A. = CA,...AAHx),dPix),ddPix),..:). (69) 

This allows, in particular, to write 

6C{x)[P]b,...b, _ ,Yi dC{x)B,...B,{P,dP,...) ^j 



in the functional differential equations (j67p and (jOSp . This completes the first step of the reduction 
process of equation ()60p . to a countable set of distributional differential equations. 

The strategy to convert these into regular differential equations now begins with eliminating all 
6 distributions, which requires to shovel derivatives over to the test function. For the zeroth order 
equation (j67|) we obtain 



= I dx\ CAU^^''{x){dl^J){x,x) - CA{x)S^^ix){d2M^,x) (71) 



j=0 s=0 



for any test function /(x, y) with compact support. Unfortunately, one cannot directly read off from 
this equation that the coefficient functions of the various derivatives of / all vanish. This is because 
the derivatives dif and 82/ of the test function are evaluated at {x, x) rather than (x, y), and thus 
are not independent of each other. Indeed, we have 5^/(x,x) = (5i^/)(x,x) + {d2fj,f){x,x), so 
that 

dx{A{x)f{x, x) + B>'{x){di^f){x, x) + C'^(x)(92^/)(x, x)] 
dx{[A{x) - 5^C^(x)]/(x,x) + [B^'ix) - C^{x)]{di^f){x,x)] . (72) 

In particular, the vanishing of the first integral for any arbitrary test function / only implies that 
A — d^C^ = and — = 0, but not that the coefficient functions A, B^^ and would vanish 
individually. This applies similarly if higher order derivatives are involved, since with 

n 
t=0 
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we can always express derivatives acting on the second entry of / by those acting on the first entry 
and total derivatives, and then read off the independent equations. Using (I73p and re-ordering 
multiple sums, the zeroth order equation (i7T|l can be brought to the form 



0= dx 



oo 

f{x,x)A{x) + J2(9r(S,...^J){x,x)B^^-^-{x)} 



(74) 



where the vanishing of the coefficient A amounts to the differential equation 



oo J 



the vanishing of the coefficient B^^ to the differential equation 

dC 



— d^-^ 



(N=0 



j=0 



oo ] 



j=0 s=0 



dC 



pA 



ai...as)7i---7j 

and the vanishing of all further coefficients i?^i^2--- ^]^g differential equations 

dC 



(N=0 
{w>2 



) 0= ^ Ll,)dil'-M^^^^ 



j=w-l 

oo j 



ggj pA 

\l32---liw)ctw-aj 



j=w — l s=w—l 



fP~^ /If 



dC 



pA 



\lil...l3u,Oll...Ols-nj)"l\---"lj 

This countable set of partial differential equations for the coefficients C and Ca is equivalent to 
the information contained in the one functional differential equations (j67p arising at order = 
in K. 

Similarly, one obtains for each order N >1 from equation (j68p first the distributional differential 
equation 

= j dx[{N+ 1)\CaB,...B^ [U^^'dl^sf - S^"'d2^f) + NN]QB,^^CB,...Br,)Md2^ 



-N\ 



oo oo y .\ n„ 

yy A ^^B,...B^ ^-s ^Ac 



s=0 



oo oo N—1 



5C 



Bi...Bi...BM Qs 



^ fioj pR 2{ai...as (i+2) as+i...aj}- 

1 i=s i=l (^(Jax-ajlr 

-^t-k dCB^...BN-i 



oo oo oo 



+(Ar-l)! V V V (-1V ^— ^ 
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where 5(3 ^y+i)/ denotes a derivative acting only on entries three to + 1 of the test function. 
The last multiple sum arises from an elimination of the partial derivatives acting on entry number 
+ 2 of the test function by way of the identity 

3 

'^{2,N+2)ai...ajf ii)'^iai...as(~^y '^(1,3,...,N+1) as+i...ajf ' C^^) 

3=0 

which renders the distributional differential equations ()75p for each N free of derivatives dis[+2f 
and thus removes ambiguities due to surface terms, so that one can now write ([75]) in the form 

/oo oo 
E E E (9id-^ ■ ■ ■ dZrh^.Mf - (^2 ^ , (77) 

s=l j=0 Part^O ) 

where the third sum is meant as the sum over partitions j = + . . . uin^i. Employing various 
multinomial distributions of higher derivatives and reordering of sums one obtains the following 
equations for N > 1. At level j = one obtains from the vanishing of the coefficient ^^'^^B the 
equation 

U,.f^f^...J = (AT + 1)\Cab,...b^S^^ - NN\Q^s,''^Cb,...b^)m 

oo N-l aft 

y^CB^^ , M^^ + (Ar-l)!y — ^^4^ 
from the vanishing of the coefficient ^'^''^^B the equation 

L=2-MUo) = {N + 1)ICas....s.U^^^^^ - NI±J J'^'-^^ di^L^M^l^^^ 

j=l ^'^{I3^\a2...a,^ 
N—1 af oo 

and from the vanishing of the coefficients ^^-^'^^B the equations 

( ) = A^'V('M ^'-^b^...Bm A|/3,) 

U>3;j=0+- +0; U ^""^^ dd^ pA 

i=s-l (/3i.../3i,-i|as---aj 

Af— 1 Qfj oo 

+(Af-i)! y - (AT - y ^. . ^^^-^^"^ 

i=l ^'^/3i.../3/ j=. ^'^i.../3,a,+i...a/ 

At level j > there are two more types of coefficients that lead to equations. The first type is 
(«>i;i="i-a+2)^^ where the a-th member of the partition ma+2 = j; and their vanishing leads to the 



29 



equations 

dC ~ 

( ^ 0- (N -1)U'^^] B^-Bg-BN 

U>l;j=0+- +j+-+0; U-l^iV s ) ^^g^j pj^^ 

and the second type (*^2i;-J=P^''t™0))iJ covers all remaining partitions of j > 2, which have at least 
two non-vanishing members, and their vanishing leads to the equations 

( N>1 N n_ {N-l)\ ^ ( gg-j-s dCB,...B^_, 

Vs>l;i=m3+-+mjv+i; m3!...m7V+i! ^ ^ ' s\{q-j-s)\ °'-+^+^-°"' QQl ^ pB^ 

q=s+j Pi...lJs+jas+j + l---0:q 

Fortunately, these equations encoding the first Poisson bracket relation [3T] considerably simplify 
upon further inspection, as we will show in the the following section, where they will also be 
supplemented by equations equivalent to the remaining second Poisson bracket relation (j32p . 



F. Construction of the local superhamiltonian 

The differential equations for the coefficients CbiB2... imply that the latter only depend on 
at most second order derivatives of the P^. For one first observes that insertion of equations 

(A^ > 1, s > 2, ma H h niN+i > 2) into the equations {N > 1, s > l,j = ^ h j + . . . > 2) 

yield 

^'^B^:-^^-;^'^ = , (78) 

first apparently restricted to > 1 , but then insertion of this result into the difference of equations 
(iV > 1, s = 2, j = 1) and (A^ > 1, s = 3, j = 1) shows that ([78]) holds in fact for ah > 0. The 
only other conclusion one may draw from the last two sets of equations of the previous section is 
that for a = 1, . . . , A^ we have the symmetry condition 



^^B^...Ba...BN _ Bn- 



Qf)2 pBa QQ2 pBM 



for all Af > 1 . (79) 



Insertion of these strong results into the remaining three sets of equations for A^ > 1 collapses the 
latter to two equations coupling coefficients of orders A^ + 1, A^ and A^ — 1, 



={N + 1)! Cab,...b^U'^''^ - N\ ^^^ki^M^I") - 2Ar! ^'^f^i-^iv 



JN -2){N (80) 
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and 



=iN + + iN- 1)! V f"--^--"^^ - 2(iV - l)!a 



^^^ Cb,...B^ ^Aa ^ ^l^^B^^Q Aa ^ ^^ dCB,...B^ 



Ma/ 



-iV7V!Q(B/*'^"CB,...B^)M- (81) 
as well as a further symmetry condition 

= ^A|7) foj. all AT > , (82) 

where the A'^ = case is provided by the equation (A^ = 0, w = 3). The only other independent 
equation is the one for {N = 0,u; = 1), coupling C to Ca, 

=2d,iCAU-^n + 2CaS-^ - 2d, (^^"'^^) - 

dpA ^ dd^^p^ ddf,,p^ ^ ' 

since the equation [N = 0, w = 0) is simply the divergence of this, and all equations {N = 0, w > 4) 
are identically satisfied. Thus only the five sets of equations (j79|) . (j80|) . ([HT]) . (j82]) and ([83|) must 
be solved for the coefficients CAi...Ai^{P-,dP,d'^P). 

But in addition to these equations, the weight-one densities CBi...Biq{P ,dP P) must also 



satisfy three additional conditions 2l| imposed by their transformation properties under changes 
of coordinates on the hypersurface S (equivalently, these follow from the Poisson bracket of the 
supermomentum and the superhamiltonian). Under an arbitrary change of coordiantes = x°^{x), 
the fields - transform as 

pP,...Pi ^ . . . {A-^f^\, (84) 

where {A~^)a = dx^ /dx^ is the inverse of the Jacobian A'p = dx"^ jdx^ of the transformation. 
Since the coefficients Cbi...Bjv [-^] ^^^^ ^ tensor densities of weight one, they transform as 

Cc,...cAP^ dP, d^P) = det(A)^g . . . A^ICb,...bAP^ dP, d'^P), (85) 

where A^ = Aif^^ . . . A^j^ denotes the transformation of the capital multi-indices. Taking the 
derivative of equation (|85]) with respect to {A~^)a,f_iu = d^x^ / {dx" dx^dx'^), noting that its right 
hand side is independent of these quantities, we obtain quite generally 

dcgP 

Q = J pa2-ai(cr ^'-'Bi...Bm /gg-v 
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This is the first invariance identity for the coefficients Ca^.^a^^ that also follows directly from 



the constraint algebra. Taking the derivative of equation ()85p with respect to [A 
d'^x''' / {dx^dx'^) and using the first invariance identity (I86p we obtain a second invariance identity: 

degP 



1=2 



J pa2...ai(fM dCBi...Bf^ _ q pai...ai dCBi...Bj^ +2/(9 p"2---Q/(A' ^^Bi...Bjsi 

QQ pa2...aip ^ QQ2 par...aj " QQ2 pa^.-ajp 

1^) pv u)a 



(87) 



The last invariance identity is obtained by taking the derivative of (j85p with respect to {A 
dx^ /dx^ which results in 

dCB^...BM I ppfe.../3r 9Cb^...Bm 

(88) 



2Q ph-0T ^^Bi...Bn 



(i) (i) 

where rii is the number of small indices contained in the capital index Bi and Bi = P\ . . . /3„ . . If 
we contract the indices we get the simpler indentity 

dcg p 

with n being the total number of lower case indices contained in all capital indices Bi to -Bat. 
Equations ([86]), ([87]) and ([88]) together with the equations ([79]), ([80]), ([HI]), ([82]) and (gSD must now 
completely determine the coefficients Cbi...Bi^- These then yield the local part of the superhamil- 
tonian, so that together with the already explicitely known non-local part and supermomentum, 
this determines the gravitational dynamics. The physical problem of finding dynamics for modified 
dispersion relations is thus reduced to the mere technical problem to solve this set of homogeneous 
linear partial differential equations. 



G. Reduction to first derivative order 



We remark that the linear partial differential equations determining the local part of the super- 
hamiltonian can in fact be reduced to linear partial differential equations for quantities that depend 
at most on the P^ and their first partial derivatives. This follows essentially from the observation 
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that the coefficients CbiB2... depend first of ah only polynomiahy on the second partial derivatives 
of the P^, and indeed at most to order diniM — 1. Since due to the fact that the coefficients 



Cbi...Bm do not depend on derivatives of higher than the second, we can ffi'st extract a further 
set of equations from ([8T]) . Writing out the total divergence of the third term, we then conclude 
that 

dCB,...B^ -0, A^>0. (90) 



ddf P^dd"^, P^ 

[per ii)u 



For transparency, we restrict the following technical discussion to the case where S is a three- 
dimensional manifold. However, the argument holds in a generalized form in any dimension. Since 
the coefficients Cbi...Bi^ are tensor denisities it can be checked that for all > 

B,...B^ Q R S T - qq2^^Pq QQ2^pRQQ2^^ps qq2^Pt ^ > 

are also components of a tensor density. According to equation (|90|) the quantities A vanish 
whenever we symmetrise over three adjacent greek indices, which also implies that the A are totally 
symmetric under the exchange of the pairs a/3, 7^, k\ and pa. Moreover, the A's are also totally 
symmetric under the exchange of Q,R,S,T. Let us now investigate all the above components. In 
three dimensions it is clear that at least three of the eight greek indices in (I9ip take the same value. 
Using all the described symmetries we can always arrange for these equal indices to appear right 
next to each other, which immediately implies that 

» 75 kA po- _ /Q„N 

^^Bi...Bn Q R S T ~ ^- K^'^J 

Put another way, in three dimesions, the coefficients Cbi...Bn can depend on the second derivatives 
of P^ only up to the third power. We may thus expand 

CB....B. = (^)Ab,...B. Q'/s'P'',a, ^^A + ^^^A^^ q' / P% (93) 

+ ^'^^B....B^Q'P'',a, + ^'^^B.-B., (94) 

where the coefficients ^*^A can depend on the P"^ and their ffist derivatives only, and only the 
highest order coefficient (^-'A must transform as a tensor density. In this way the dependence of 
the coefficients Cbi...Bi^ on the second derivatives of P^ can be completely eliminated from our 
differential equations. If S is of higher dimension we simply have to add more derivatives in (|9ip . 
Thus, in general, the coefficients Cbi...Bn depend polynomially on the second derivatives of P^ at 
most to order dimS. The coefficients now have to be determined from the remaining equations. 
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H. Example: Canonical dynamics of second degree dispersion relations 

We now illustrate how to solve the linear partial differential equations we identified in sec- 
tion IIIIFI in order to obtain gravitational dynamics, for the simplest case deg-P = 2. On a 
four-dimensional manifold, this directly yields Einstein-Hilbert gravitational dynamics with un- 
determined gravitational and cosmological constants (which appear as integration constants and 
must be fixed by experiment) as was first shown in [8] a long time ago. The point here is of course 
that we have the relevant equations for any admissible dispersion relation, not only those of second 
degree, and only wish to illustrate that one can indeed proceed from these equations without fur- 
ther assumptions in order to obtain the gravitational dynamics of the specific spacetime geometry 
at hand. In particular, due to our foregoing comprehensive analysis that extracted all information 
from the constraint algebra, we do not need to draw on any results beyond our equations. 

In the case of a second rank tensor field P, which we consider here, all capital indices con- 
tain symmetric pairs of lower case greek indices running from 1 to 3. First we observe that the 
coefficients M^^ and Qb^^ vanish since the non-local part of the Hamiltonian is equal to zero. 
Moreover, the coefficients U^'^^ and S^^ reduce to 

Thus equation (|80l) for iV = 2 reads 

= Cab.b^U"^^^ , (96) 

which can be directly solved yielding CaBiB2 = 0. Inserting this result into ([80]) . starting with 
= 4 and iterating on all even A^, we find that all coefficients Cbi...Bm with an odd number 
of capital indices greater or equal to three already vanish. For calculational convenience only, we 
perform a change of variables from P"^ to with P'^'^g-yp = 5^ and substitute the first and second 
partial derivatives of P'^^ by those of ga/s accordingly. After this change of variables, equations 
(1861) and ^ become 

^5^^-^" = and _^!%-^i^ = 0. (97) 

Using a similar argument as in the previous section, we can now show that we even have 

/^%^ = 0, (98) 

because, in three dimensions, either one of the indices 1,2,3 appears at least three times, so that 
all components of these tensor densities of weight one vanish according to the above symmetry 
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conditions. Thus all remaining coefficients can depend at most linearly on the second derivatives 
of gai3- This has the direct consequence that according to equation (ISOh with = 1 the coefficient 
Cab cannot contain second derivatives of gap- But then equation (f80|l implies Cbi...Bm — for all 
even > 4. Hence, it remains to determine Cab, Ca and C to find the gravitational dynamics. 
We start with the discussion of the coefficient C. We already know that it has to be linear in the 
second derivatives of gai3 so that 

C = Co{g,dg)+Af^'gap,^s, (99) 

where Ai is a tensor density of weight one and contains no second or higher partial derivatives of 
gap- However, equation (f87|) implies that Ai, being a tensor density, can not even depend on the 
first partial derivatives of y^^. Thus equation ([89|) for A^ = can be rewritten into 

C = ^R^fs-ysAf^' + Ao{g), (100) 

where Ra/s-yS is the Riemann-Christoffel tensor of the metric gap and Aq a tensor density of weight 
one that is solely constructed from ga/s- In three dimensions, the Riemann tensor can of course 
be expressed in terms of the Ricci tensor Rai3 = P'^'^R^aSp- Now ^/^^^^etg R, where R = RapP"^ 
denotes the Ricci scalar, is the only weight-one scalar density linear in the second derivatives of gap 
that one may construct from the Ricci tensor and gap-, and the only scalar density of weight one 
one can construct from gap alone is \/— det g. The minus sign under the square root accounts for 
the fact that with our normalisation condition the metric on the hypersurface is negative definite. 
Thus, we finally arrive at 

C = -(2K)~V-det5(i?-2A) (101) 

with two real integration constants k and A. It is then simple to determine Cab from equation 
dlOD for A^ = 1 and we find 

Cap^lu = (16k)"^ yj- det g [ga^iOpu + 9p^,9au - '^.gapQtiA ■ (102) 

Finally, we calculate the coefficient Ca, which depends at most on the second derivatives of gap 
and is at most linear in those. Equation ([83]) reduces to 

= P'^^P^^VaCp^, (103) 

where we use the the torsion-free covariant derivative Vq compatible with gap only for notational 
convenience. Using equations (|86p - ()88p and following a similar argument (2^ as for the coefficient 
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C immediately yields 

C^p = /3iV-det5 {Ro,p - l/2gapR) + h^J-deigg^p (104) 

for some constants /3i and /32- The remaining equations (I79p with N = 2 and ()8ip with N = 2 
involving Cq/? are then identically satisfied. We note that the coefficient Cap can be written as the 
functional gradient 6S/6P°'^ of the scalar density S = det gR—2fi2\J — det^i, and finally make 

the transition from the full Lagrangian (j65|) to the superhamiltonian by means of the Legendre 
transformation (j5ip - ()53p . For the canonical momenta T^ap^ one then has 

v^„, = ^ = 2C.,,,i^-^ + ^. (105) 

However, the canonical momenta (j36p are only determined up an additive functional derivative of 
some scalar density of weight one with respect to P°^^ . One can thus drop the second term on the 
right hand side of (jlOSp by redefining Ti^p T^ap ~ ^S/5P'^I^ without changing the dynamics of the 
theory. Then the superhamiltonian reads 

H = C'^^^T^apT^-^s - 2C''P^^CapT^^s + (2K)-V-detg(ii - 2A) (106) 

with 

(jap^S _ ^ (^pa"t pPS _|_ pl3^ paS _ paP p^S-^ (107) 

\/— det g 



due to 



The second term in the superhamiltonian can be shown to be dynamically irrelevant 
the special form of the coefficient Cap- With the superhamiltonian (|106p and the supermomentum 
Da = 2p^'^'V pffa-y from (j44p . we have finally found (as did for a construction that only works 
for degP = 2), the gravitational dynamics in the case of a three-dimensional hypersurface S for a 
hyperbolic polynomial of degree two, also known as general relativity. The task to find canonical 
dynamics for dispersion relations beyond second degree is now of course to find solutions of our 
equations for degP > 2, which appears a much harder task. But this precisely what it takes if one 
wishes to consider modified dispersion relations in earnest. 
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IV. DYNAMICS OF TENSORIAL SPACETIMES 

We finally address the master problem of deriving the equations determining the gravitational 
dynamics of a fundamental geometric tensor field G, under the assumption that the latter gives rise 
to a hyperbolic, time-orientable and energy-distinguishing tensor field P by virtue of specific matter 
field equations. This gravitational theory for G is more fundamental than the phenomenological 
gravity dynamics derived for P in the previous chapter, since fields couple directly to G, and so 
do point particles via P constructed from G. But this greater generality comes at the price of a 
less sweeping construction scheme. While the always totally symmetric, even rank tensor fields P 
can be treated in precisely the same way for any rank, the fundamental geometric tensors G come 
in all possible ranks and symmetries (as long as one can couple matter fields to them), and the 
construction of their geometric phase space must proceed in fashion of a case-by-case analysis. But 
apart from these technical details, the overall construction is as simple and inevitable as in the 
previous chapter, and one obtains also a system of homogeneous linear partial differential equations 
whose solutions determine the gravitational dynamics of the geometry (M,G). 

A. Construction of tensorial spacetime geometries and their dynamics 

The construction of gravitational dynamics for a fundamental geometric tensor field G proceeds 
logically exactly along the same lines as that for dispersion relations. The only relevant difference 
consists in the choice of canonical variables for the dynamics, and all the technical modifications this 
entails. To separate the essential steps from their technical details, we therefore quickly prescribe 
the general recipe one has to follow to make any candidate geometry (M, G) into a spacetime 
structure and to derive the equations determining their canonical dynamics. We will then see this 
recipe in action in the next two sections. 

1. Decide on a tensor field G of arbitrary valence as the geometry on a smooth manifold M. 

More generally, one may also choose a collection of tensor fields G = (Gi,G2, ■ ■ ■), each 
possibly of different valence, to provide the geometry. This would be the case, for instance, 
if one aimed at studying a bosonic string background featuring a metric g, a two-form field 
B and a scalar cf), say, using the philosophy of this paper. 

2. Decide on matter dynamics to define the causal structure impressed on the geometry. 
These matter equations may well be of phenomenological nature, as were the Maxwell equa- 
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tions before 1905. This is where the theory gets its vital injection from reahstic physics. 

3. Calculate the totally symmetric covariant tensor field Pq associated with the linear(ized) 
matter field equations in terms of the geometry G. 

This is straightforward in principle, but may in practice require to first remove gauge am- 
biguities. If several matter field equations are present in the theory, one needs to consider 
their entirety to calculate P. 

4. Restrict attention to geometries (M, G) for which Pq is hyperbolic, time-orientable and 
energy-distinguishing. 

Only these deserve to be called spacetimes. In order to get an overview over which algebraic 
classes of the geometry G present spacetimes, it is often useful to figure out the algebraic 
classification of the geometric tensor G under GL(dimM) transformations and associated 
normal forms. 

5. Construct the configuration variables describing the spatial geometry on an accessible initial 
data hypersurface by normal and tangential projections, eliminate the degrees of freedom 
fixed by Pq = 1 and Pq = 0, and associate canonically conjugate momenta to all remaining 
degrees of freedom. 

6. Construct the supermomentum and superhamiltonian exactly along the same lines as done 
in section IIIH but with the spatial point particle geometry replaced by that for fields, as we 
will illustrate for area metric spacetimes in four dimensions in the following two sections. 

The precise form of the coefficients will depend heavily on the geometry G chosen. But once 
a concrete geometric tensor G is chosen, and its independent degrees of freedom have been 
identified, the calculation goes through also in this case without complications. 

7. Solve the resulting system of linear partial differential equations to determine the local part 
of the superhamiltonian. 

How difficult this is now very much varies with the geometry G that has been chosen. 

For the simple case of metric geometry carrying Maxwell theory, execution of this programme leads 
to the condition that the metric must have Lorentzian signature, and the system of homogeneous 
linear partial differential equations has a unique family of solutions, giving rise to the standard 
Einstein-Hilbert gravitational action with undetermined gravitational and cosmological constants 
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(which appear as constants of integration). Essentiahy, this has been shown a long time ago [s], 
and is of course recovered as a very special case of our general construction. 

Any other tensorial geometry requires a separate case-by-case analysis for virtually all of the 
above steps. We therefore choose to illustrate the procedure for area metric geometry, which 
accompanied us throughout this paper as a particularly interesting example for the workings of 
our general theory. 



B. Phase space for area metric geometry seen by electromagnetic fields 

To illustrate the procedure of finding canonical dynamics directly for a fundamental geometric 
tensor field G underlying a chosen field theory, we will concentrate, for definiteness, on the partic- 
ular example of a four-dimensional area metric geometry coupled to electromagnetric fields. We 
assume that the inverse area metric G"'^'^'^ is everywhere non-cyclic such that with ([4]) the totally 
symmetric geometric tensor p^'^'^ takes the form 

24 

p9^bcd ^ Qmnr{aQb\ps\cQd)qtu_ 

^ {G'l^kijklY 

Using the complete spacetime covector frame {n,e^,e^,e^} constructed from P along a hypersur- 
faces X{Ti) given in terms of the embedding map X{y), we then define the functionals 

= G(n(y),6-(y),n(y),e''), (109) 
G%(y)[X] = ^u;G^^sG{niy),e^{y),e^y),e'{y)), (110) 

G^piy)[x] = \u;Gaj5i^Gp^.M^ny),e\y),e''{y),e''iy)) ■ (m) 

where we used the volume form uJcap-y = (— det G°^)~^''^eQ/3^ induced by the symmetric hyper- 
surface tensor field G°^^ to construct the hypersurface tensor fields G"^ and G^p from the other 
possible projections of the inverse area metric onto the hypersurface X(Ti). Note that the index 
positions really distinguish unrelated tensor fields, which together encode the degrees of freedom of 
the inverse area metric on the hypersurface. From the normalisation conditions PG{n,n,n,n) = 1 
and PG{n,n,n,€") = it follows that G"^ can be assumed to be trace- free and symmetric with 
respect to G°'^ . 

The phase space of a four-dimensional area metric spacetime is then spanned by tensor fields 

G"^ (112) 
on the three-dimensional manifold S and their respective canonical momenta 

n„^, n/, n-^ (113) 
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which are taken to be tensor densities of weight one. Adjoining canonical momenta again is then 
equivalent to adopting a Poisson bracket 



{E,F] = j dy 



6E 5F 6F 6E 



5G^ SUa 6G^ SUa 



(114) 



on the space of functionals of the canonical variables and fl^, where the capital index A 
collectively denotes the different greek indices with their respective positions: = ("^,°^,a/3)- 

Again we will look for dynamics in terms of a Hamiltonian as it appears in (j27p that evolves the 
phase space variables (G"^, 11^) with an evolution parameter t, such that the embedding of the data 
at time t by virtue of a foliation Xt : T, ^ M produce an inverse area metric G"^^*^ on M whose 
dispersion relation is hyperbolic, time-orientable and energy distinguishing. The 21 components 
of the inverse area metric can then be reconstructed from the 17 independent components of 
the symmetric tensor field G"^, the trace- free hypersurface tensor field G"^ which is symmetric 
with respect to G°^^ , the symmetric hypersurface tensor field Gap the spacetime vector frame 
{Tt,eta} by 

G'^'^'^iXtiy)] = AG^'Tte%Tl'ef, + G,,.(a;^i)''"^(a;^i)-^^e?„e,%e?,ef, 

+2(Gf ^ + 6^ ^){u-'r'Tl''e%et,efs. (115) 

The conceptual steps in the construction of the supermomentum T> and the superhamiltonian H 
on the phase space given by (jll2p and (jllSp are precisely the same as for the pure point particle 
geometry in section IIIII We will quickly go through these steps in the next section. 



C. Canonical dynamics for area metric spacetime 

We already saw in the previous section that the canonical phase space in the case of a four- 
dimensional area metric spacetime consists of the tensor fields G"'^ (symmetric), G"^ (trace- free, 
symmetric with respect to G"^) and Gap (symmetric) as well as their conjugate momenta flap, 
fla'^ and fl"'' with the same respective algebraic properties. The superhamiltonian 7i{N) and the 
supermomentum T)[N°^da) satisfy the Poisson algebra relations ()3ip . ()32p and ()33p . but now the 
symbol Pq^ on the right hand side of (j3ip is not a canonical variable itself, but the particular 
phase space function 

p^^ = 1 [g'^^G^'^G^s - G'^^G^^G^s - 2G°^G^ ^G^ ^ + SG^^G" ^G^ ^) . (116) 
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The construction of the supermomentum follows the same steps as in the case of the point particle 
geometry. One readily finds 



dyN^y) (a^G"^)n„;3 + 29«(G"^n^^) + (a^Ga/3)n°^-29«(fl°^G^^) (117) 



(118) 



The non-local part T^non-iocai of the superhamiltonian can be found using the first of the necessary 
relations (f29]) . One calculates 



n{N) 
njN) 
n{N) 



N{z)[...]+d,N{z) [-2(u;^i)^^("G^V 



N{z)[...]+d^N{z) 



N{z)[...]+d^N{z) 



GpaT Lr 



iu;^y"'^G^p 



n"^(z) 

which may be integrated to yield the non-local part of the superhamiltonian 



n 



non-local 



2d^ 



(y) + 39-, 



-G.ri,G\)Ppn'^^ 



(y) 



(119) 
(120) 
(121) 

(122) 
(123) 



It remains to evaluate the Poisson bracket of two superhamiltonians to determine its local part. 
The analysis proceeds along precisely the same lines as in the case of the pure point particle 
geometry. We perform a Legendre transformation of the local part Hiocai of the superhamiltonian 
with respect to the momenta IIa- The resulting Lagrangian L[G'^]{K'^) then satisfies the linear 
functional differential equation 
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The coefficients M^/^ and can be read off from the functional derivatives of the non-local 

part of the superhamiltonian with respect to the canonical variables. The coefficients U^'^^ read 

jjalB-yS _ _Qpi'r\iaQl3)\S) ^ (125) 
jja^-yS ^ _3pa(7(A'5)^ ^ SP-^^^^^^ ^ , (126) 

(127) 

The coefficients S^'^ can be calculated from 

= O/^q^^^^'m^i^) - Q^^t'^'a^^M^i^] - dfsU'^'^"' - w^^dpG^ - w^^, (128) 
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with = 2P^(°95G^)^ ^ = P^'^dsG^ ^ - P^^d^G'^s and V^/^ ^ = -2P^^d^^^Gp)5. Expand- 

ing the Lagrangian L[G'^]{K^) into a power series in the velocities K^^ 

oo 

L{x)[G\K^) = Y,C{x)[G%,...B^K^-...K^% (129) 

i=0 

one derives exactly the same equations ([79]) - ([83]) for the coefficients Cbi...Bn as in the point particle 
case. The coefficients Cbx...Bn ^.re again tensor densities of weight one and, as a result of the 
algebra equations, depend at most on the second partial derivatives of the fields G^. Since the 
hypersurface S is of dimension three, it again follows that the coefficients depend on the second 
partial derivatives only up to the third power. 

The invariance equations following from the transformation properties of the weight-one den- 
sities Gbi...Bm (o^i fully equivalently, from the Poisson bracket of the supermomentum with the 
superhamiltonian) can be derived in straightforward fashion. The first invariance identity takes 
the form 

^ ^^.jg dCB, ^jg dCB, _ ^ OCb, B, _ OGb, B, 

ddl^^Gpp 'ddj^f^ ''ddl^G,^, 

The second invariance equation reads 

g 2(^m(" 9Gbi...Bi^ _I_ QfJ'ja ^Gbi...Bjsi _ Q Qt^'^ dC Bx...Bis[ 

Q{a dGBi...BM _ Qti dGBi...BM , 2d ^Gbi...Bm _ 2Q QfJ. ^Cbl^Bm 

'dd,)G\ 'dd^^G%^ 'dd}^G\ ^ 'dd\G%^ 

o dGBi...BM or< ^Gbi...Bn .0 a dCB^...BN r. ^ dCB^...BN ^-.o-.^ 

ddlpG^u ddi^^Gp)p 5^'(a^/3)M ddlpG^,, 

The last invariance identity is even more complicated and we only display its contracted form: 

Ci^^^r — nr<l^iy ^^Bi...BN dCBi...BN 

- (3 + n)CB,...B^ - 2GP - 2G,.^^ 

Q^,. dCB, B. _ ^p 9Cb, _ ^ OGb, B. 

_ 2dl^G^^ ,^^^1^ - AdlG'^''^^^^, (132) 
""dd^G^, dd%G^,' 

where n is the difference of the total number of subscript indices and the total number of superscript 
indices in the coefficients Gbi...Bn- 



42 



V. CONCLUSIONS 

In this paper, we addressed the question of how to construct canonical gravitational dynamics 
for spacetime geometries beyond the Lorentzian manifolds featuring in Einstein's general relativity. 

The first step consisted in an analysis of what kind of tensor fields G on a smooth manifold 
M can serve as a spacetime geometry in the first place, dependent on the presence of specific 
matter field dynamics. Indeed, the geometry must be such that all matter field equations are 
predictive, interpretable and quantizable. These conditions on the matter field dynamics impose 
three corresponding algebraic conditions on a totally symmetric tensor field P, which is defined in 
terms of the fundamental geometric tensor field G and whose precise form arises from the matter 
field dynamics: P needs to be hyperbolic, time-orientable and energy-distinguishing, as reviewed 
in section III A[ So in order to start the whole machinery presented here, we first need to know 
which matter couples in which way to the tensorial geometry. We do not see this as a weakness of 
the formalism, but rather as an insight; it was Maxwell theory that justified Einstein to promote 
Lorentzian manifolds to the status of a spacetime geometry, and experimental observation of any 
matter that does not mimick the structure of Maxwell theory (non-half-integer spin or superluminal 
matter, for instance) will force us to choose another tensorial geometry. But certainly one that is 
hyperbolic, time-orientable and energy-distinguishing. Fortunately, with the results of chapter |TIl 
we have all these geometries under excellent technical control. 

Directly from these purely kinematical insights, one can calculate the deformation algebra of 
hypersurfaces in any hyperbolic, time-orientable and energy-distinguishing geometry. This is the 
algebra of linear operators that describe how the geometry induced on a hypersurface changes when 
the hypersurface is deformed in normal and tangential directions. And gravitational dynamics is 
precisely about understanding these changes in the geometry on initial data surfaces, as has been 
clarified in seminal work of Hojmann, Kuchar and Teitelboim for the special case of Lorentzian man- 
ifolds, building on the canoncial formalism introduced by Arnowitt, Deser and Misner. The most 
important result of chapter |II1 from a practical point of view, is therefore that one can calculate the 
deformation algebra of hypersurfaces in any hyperbolic, time-orientable and energy-distinguishing 
tensorial geometry. This is by no means trivial, since the existence and uniqueness of the way to 
associate normal directions along a hypersurface with its canonical normal co-directions by means 
of a Legendre map requires all three algebraic properties: the hyperbolicity, time-orientability and 
energy-distinguishing property. Despite appearances, this also applies to metric geometry (where 
any one of these conditions is equivalent to the requirement of a Lorentzian signature), since al- 
though, purely formally, one can still construct normal directions from normal co-directions for 
other signatures, they lose their physical meaning. 
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Canonical gravitational dynamics for the spatial geometry that are ultimately invariant under 
spacetime diffeomorphisms must be given by a pure constraint Hamiltonian (which is of course a 
functional of the geometric degrees of freedom and associated conjugate momenta on an initial data 
hypersurface) composed of two separate first class constraints — corresponding to spatial diffeomor- 
phism invariance within the hypersurface on the one hand, and invariance under diffeomorphisms 
away from the hypersurface on the other hand. The role of the deformation algebra, in this geo- 
metrodynamic language, is that these constraints must satisfy canonical Poisson bracket relations 
of the same form as the commutator algebra of the normal and tangential deformation operators 
on functionals of the hypersurface embedding map. The task is thus to determine the constraint 
functionals from this Poisson algebra. 

This would be a mere representation theory problem if the Poisson algebra were a Lie algebra. 
But only two of the three bracket relations feature structure constants, and their impact on the 
form of the constraint functionals is thus readily established. The third Poisson bracket relation, 
however, features a structure function that captures the impact of the particular hyperbolic, time- 
orientable and energy-distinguishing tensorial geometry to which one wishes to give dynamics. 
Determining the implications of this bracket amounts, at first sight, to the truly daunting task 
of solving a system of non- linear functional-differential equations. The better part of chapters IIIII 
and IIVI is thus devoted to reducing this to the equivalent, and principally manageable, problem 
of solving a system of homogeneous linear partial differential equations. And this set of equations 
already contains, by construction, all possible classical gravitational dynamics one can give to 
a tensorial spacetime geometry that can carry predictive, interpretable and quantizable matter 
fields. In the philosophy of this paper, the physical problem of finding diffeomorphism-invariant 
gravity theories alternative to Einstein's general relativity is shown to be equivalent to the mere 
mathematical task of solving these linear partial differential equations. 

This casts important physical questions into precise mathematical form. The question whether 
there are any alternatives to general relativity turns into the problem of existence of solutions; 
the question whether there is a choice between various dynamics for a given tensorial spacetime 
geometry translates into the question of their uniqueness; and finally, the actual construction of 
all concrete gravitational dynamics amounts to nothing more, but also nothing less, than explicitly 
finding the exact solutions of these linear partial differential equations. 

The difference between the treatments in chapter IIIII and IIVI is that only in the latter do we 
construct dynamics for the fundamental tensorial spacetime geometry G to which also fields can 
couple, while in the former we give dynamics only to the totally symmetric tensor field P seen 
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by point particles. In the special case of Lorentzian geometry, the two points of view accidentally 
coincide, since the tensor field P encodes precisely the same degrees of freedom as the fundamental 
Lorentzian metric g to which fields can couple. The key result of chapter IIIII is that there is at 
most a one-integer-family of essentially different gravity theories that differ in their prediction of 
particle trajectories, and we wrote down the complete set of equations for all these theories. So if 
one is interested in the motion of massive and massless point matter only, one can ignore which 
particular fundamental geometric structure underlies a hyperbolic, time-orientable and energy- 
distinguishing dispersion relation, and compare observational data with these phenomenological 
theories. In contrast, if one wishes to consider a full gravitational theory to which both point 
particles and fields can couple, one needs to construct these along the lines laid out in chapter 
IIVI The resulting theories are more fundamental, but this comes at the price that the equations 
yielding their dynamics depend more heavily on the specific type of tensor field G and require the 
explicit specification of the predictive, interpretable and quantizable matter coupling to it. Our 
derivation of the relevant equations for four-dimensional area metric geometry carrying general 
linear electrodynamics at the end of chapter IIVI however, shows that also this more fundamental 
programme can be executed. 

The main open question is how to find solutions to the system of homogeneous linear partial 
differential equations in either case. But this will well be worth the effort, since solving these 
equations immediately allows to answer a string of pertinent physical questions in gravity theory. 
Four questions of high relevance for a number of current research programmes are how to 

(i) settle the issue of which modified dispersion relations are admissible, and how they are 
determined dynamically, in order to conduct a focused search for experimental signatures. 

(ii) provide canonical dynamics to one's favourite candidate of a tensorial geometry without 



(iii) free the evaluation of observational raw data from the confines of Lorentzian geometry and 
Einstein dynamics, in favour of the wider framework that includes all hyperbolic, time- 
orientable and energy-distinguishing tensorial geometries, which one is led to consider in the 
light of matter dynamics that would qualify as non-causal in Lorentzian spacetime. 

(iv) reveal all possible classical limits of quantum gravity theories where the fundamental geo- 



Our future progress on these questions thus hinges on solving the equations derived in this work. 



further assumptions, starting from matter dynamics coupling to this geometry 




metric structure can be expressed in terms of tensor fields 
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Since not even all degrees of freedom in the area metric G enter the field P, see e.g. [13j. 
This dependence from specific matter equations is not a weakness of the approach, but rather an insight: 
there are no viable or non-viable spacetime structures as such, but they prove their validity by enabling 
specific field equations to be predictive. The causality of Lorentzian manifolds is in fact the causality 
of Maxwell theory, and thus of all the other matter dynamics modelled in the image of Maxwell theory. 
This choice is the general tensorial analogue to a choice of 'mainly minus' signature (H — • • • ^) in the 
special case of metric geometry. 

If there are such superluminal particles, they can radiate off massless particles until they are intralu- 
minal, see 

Employing the techniques for general tensorial spacetime geometries instead of those valid only for the 
special case of Lorentzian manifolds, this section follows closely the philosophy and calculation of [sl- 
The minus sign in the algebra equation involving two normal deformation operators is due to our 
normalisation condition Pin) — 1. In most standard texts on the hypersurface deformation algebra of 
Lorentzian manifolds the normalisation is chosen to be P{n) = — 1, which results in a plus sign in the 
first algebra equation. 

The change of the overall sign on the right hand side of the Poisson algebra equations in comparison 
to the hypersurface deformation algebra is due to the action from the left of the deformation operators 
on hypersurface functionals as oppposed to the insertion of the supermomentum and superhamiltonian 
in the right slot of the Poisson bracket when acting on a phase space functional. 

In the standard treatments of canonical general relativity one chooses the metric gap and its con- 
jugate momentum Tr^andard ^^"^ phase space variables. It can be checked that changing the con- 
figuration variables from the metric to its inverse results in the canonically conjugate momentum 
TTjj^ = ~5'"75^*''^standard' ^hich cxplains thc different sign in the supermomentum. 
This certainly presents an entirely new angle on — and more importantly: physically well-founded ap- 
proach to — the construction of gravitational dynamics for a non-symmetric metric or for string ge- 
ometries determined, for instance, by a metric a two-form potential B and dilaton from a purely 
canonical point of view. 

For classical theories based on a connection formulation, this may require a reformulation in terms of 
purely tensorial objects whenever possible. Possibly some pure connection theories therefore escape the 
formalism presented here, but in this case we feel a careful evaluation of their kinematical apparatus, in 
particular the definition of observer frames, should yield the relevant information. This then certainly 
requires a case-by-case analysis. 



